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Abstract 

We determine the growth rate of linear instabilities resulting from long- wavelength transverse per- 
turbations applied to periodic nonlinear wave solutions to the Schamel-Korteweg-de Vries-Zakharov- 
Kuznetsov (SKdVZK) equation which governs weakly nonlinear waves in a strongly magnetized cold- 
ion plasma whose electron distribution is given by two Maxwellians at slightly different temperatures. 
To obtain the growth rate it is necessary to evaluate non-trivial integrals whose number is kept to 
minimum by using recursion relations. It is shown that a key instance of one such relation cannot be 
used for classes of solution whose minimum value is zero, and an additional integral must be evalu- 
ated explicitly instead. The SKdVZK equation contains two nonlinear terms whose ratio b increases 
as the electron distribution becomes increasingly flat-topped. As b and hence the deviation from 
electron isothermality increases, it is found that for cnoidal wave solutions that travel faster than 
long- wavelength linear waves, there is a more pronounced variation of the growth rate with the angle 
6 at which the perturbation is applied. Solutions whose minimum value is zero and travel slower than 
long-wavelength linear waves are found, at first order, to be stable to perpendicular perturbations 
and have a relatively narrow range of 6 for which the first-order growth rate is not zero. 

1 Introduction 



In Part 1 ( Allen et alj|2007 ) we considered solitary wave solutions of a modified version of the Zakharov- 



Kuznetsov (ZK) equation which, in a frame moving at speed V above the speed of long-wavelength linear 
waves, takes the form 

u t + (u + bu 1/2 - V)u x + V 2 u x = (1.1) 

where the subscripts denote derivatives. We referred to this equation as the Schamel-Korteweg-de Vries- 
Zakharov-Kuznetsov (SKdVZK) equation as it contains both the quadratic nonlinearity of the KdV 
equation and the half-order nonlinearity of the Schamel equation. The equation governs weakly non- 
linear ion-acoustic waves in a plasma permeated by a strong uniform magnetic field in the ir-direction. 
The plasma contains cold ions and two populations of hot electrons, one free and the other trapped by 
the wave potential, whose effective temperatures differ slightly. In u is proportional to the elec- 

trostatic potential, and b = (1 — T e f/T e %) / 'y^r where T e f and T e t are the effective temperatures of the 
free and trapped electrons, respectively. As b increases, the electron distribution becomes less peaked. 
A flat-topped distrib ution is in acco rdance with numerical simulations and experimental observations of 
collisionless plasmas (|Schamellll973l ). For further background to the physical basis and applicability of 
the SKdVZK and related equations, reference should be made to Part 1. 

The existence of planar solitary wave solutions to the SKdVZK equation and their stability to trans- 
verse perturbations was addressed in Part 1. In this paper we turn to the study of planar cnoidal wave 
solutions to the equation. In Sec. 2 we show that a number of families of cnoidal wave solutions to the 
one-dimensional form of (|1.1[) exist, but not all can be expressed in closed form. Linear stability analysis 
of periodic solutions of the SKdVZK equation with respect to transverse perturbations is carried out in 
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Sec. 3. Such an analysis has been carried out on cnoidal wave solutions of the ZK and SZK equations 
which contain single quadratic and half-order nonlinearities, respectively (llnfeldl l 19851 : iMunro and Parkesl 
1999). However, as far as we are aware, such a calculation has not been performed before on an equation 
containing two nonlinear terms. The stability analysis leads to a nonlinear dispersion relation in the form 
of a cubic equation whose coefficients are finite-part integrals involving the unperturbed solution and its 
derivative. As the solutions contain elliptic functions the integrals are non-trivial. Recursion relations 
between the integrals are derived in order that only the simplest finite integrals need to be evaluated 
directly. For some types of solution it is shown that one instance of a recursion relation cannot be used 
and an extra integral must be found directly. In Sec. 4 we examine how the first-order coefficient of the 
growth rate found from the nonlinear dispersion relation depends on the type of cnoidal wave, the angle 
at which the perturbation is applied, and b. Our conclusions are presented in the final section. 



2 Cnoidal wave solutions 



To look for planar cnoidal wave solutions of permanent form travelling at speed V above the long- 
wavelength linear wave speed we drop the t, y and z dependence in Integrating once then gives 



+ Vu - -bu 3/2 



Uxx — —,.-. ., 

and multiplying by 2u x and integrating once more yields 



(2.1) 



= C + Ci 



Vu 2 - —bu b l 2 
15 



3" 



(2.2) 



where Co and C are integration constants. Although from phase plane analysis it is clear that a number 
of families of periodic nonlinear waves exist, only when Co = can closed- form solutions be obtained in 
general. Sketches of (|2.2p for the various cases leading to periodic solutions when Co = are shown in 

Fig.ru 






Figure 1: (u,u 2 ) -sketches of (|2 .2[) with Co — showing existence of families of periodic wave solutions: 
(a) C < 0, V > or b < (or both); (b) C > 0; (c) C > for some V < and b < 0. 

After introducing the variable r = s/u, (12. 2| with Co = reduces to 



Ar 2 x = g(r) = h(r) + C, 



Mr) = Vr 2 - —r 3 - V. 
W 15 3 



(2.3) 



Possible forms of g(r) for various V, b and C are sketched in Fig.O The it 1 / 2 term that appears in (II. ip 
must be interpreted as the positive square root and as a result we must restrict solutions of (12. 3p to r 0. 
In view of this, at first sight it would appear that nonlinear wave solutions to (|2.3p that cross the line 
r = with positive r 2 would have to be discarded. However, since u 2 = r 2 g(r), the r J? part of such a 
solution forms a complete closed loop that touches the origin in the (u^O, u^-plane, as can be seen to 
occur in Fig. |TJb) and (c), and hence corresponds to a nonlinear wave solution with a minimum value of 
zero. Although for such solutions in the (r, r x )-plane r x jumps from a negative to an equal and opposite 
positive value at r — 0, it is easily shown that u and its derivatives are continuous there. The jump in the 
solution in the (r, r^-plane means that the solutions are most simply expressed as functions extended by 
periodicity. Such solutions have already categorized f or the Schamel eq uation (which contains the single 
half-order nonlinearity) in lO'Keir and Parkesl ( 1997 ) . Schamell (|l972h described them for the current 
equation but in the following they are presented in a more unified form. 
The quartic g(r) will always have a stationary point at r = and also at 



'3&V 3V 3& 



(2.4) 
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Figure 2: (r, r 2 )-sketches of ^ for (a) V > 0, b > (b) V > 0, b < (c) -66 2 /25 < V" < 0, 6 > 
(d) — 66 2 /25 < V < 0, b < 0. Solid, dashed and dotted curves give rise to solitary wave, periodic nonlinear 
wave, and constant (linear limit) solutions, respectively. 



provided that V > — 6fo 2 /25. From the sketches of g(r) in Fig. [2l it is apparent that if r + is real and 
positive, nonlinear wave solutions will only occur if C > —h(r+) and the linear wave limit corresponds 
to C = -h(r+). 

If g(r) has four real roots, n < r 2 < r$ < then solving (|2.3|) yields the cnoidal wave solution 

u(x) = [r(x)\ = —— — --— (2.5) 

\ 1 + psn z (r)(x — xo)\m) ) 

where 



r<i-r 3 r 2 -n (u - r 2 )(r 3 - n) 

P= , m= p, T}=\ — , 

r3 — n ri — r 2 V 48 

and Xo is an arbitrary phase. We will refer to this class of solution as being of type I. Note that when 
V > 0, in the soliton limit c=0we have r$ = r 2 = and (|2.5p then reduces to the conventional solitary 
wave solution given in Part 1. 

When b < 0, it can be seen from Fig. [21(b) and (d) that there are periodic wave solutions corresponding 
to a g(r) with just two real roots for —h{r+) < C < — h{r~) when V > 0, and for — h(r+) < C < when 
— 66 2 /25 < V < 0. If the two real roots are r u and n, with r u > n, and the remaining two complex 
conjugate roots are a ± i/3, then solving (|2.3p in this case gives what we call the type II solution, 

, lM _ r_ M1 2 _ f Cjn + Br *) - ( Ar ^ ~ Br *) cn (^ - x o)\™) A 2 (o ^ 

u(x) - [r(x)l (A + B) - (A — B) cn(fj(x - x )\m) ) ' [ ^ 

where 
and 



A = V(r u -a) 2 +/3 2 , 5 = ^(n - a) 2 + /3 2 , 



(r„ - n) 2 - U - BY AB 

m = ; , 77 = \ / . 

4AB ' ' V 12 

We now turn our attention to the solutions written as periodically extended functions. For V > 
these occur when c > 0. As in these cases there are only two real roots, these solutions take a similar 
form to (|2.6|) . However, due to the jump in the (r, r^-plane they must be written in the form 

/ (An + Br u )- (An - Br u ) cn(fjx( X /fj)\rh) \ 2 _ = _ x f An + Br u \ 
U(X) \ (A + B)-(A-B)cn(fjx(x/v)\fh) ) ' X Cn - Br J l ' J 

where 

x(p) = (x — xq + p mod 2p) — p. 
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These solutions have a period of 2x/?7 and for a given value of V and b have a larger amplitude than the 
solitary wave. We call these type II po solutions. When 6 = they reduce to an ordinary KdV equation 
cnoidal wave solution with a minimum value of zero. 

From Fig. [2lc) and (d) it is clear that when V < 0, smaller amplitude solutions that touch u — are 
also possible. When there are two real roots, (|2.7[) still applies, while if b > and there are four real 
roots, what we will refer to as the type I po solution results. It is similar to (|2.5[) but must be written as 

. _ / r 4 + ripsn 2 (f]x(x/v)\m) 
V 1 + psn 2 (r]x(x/v)\ m ) 

When there are four real roots and b < 0, the solution which has a minimum value of zero is the same as 
the above after making the interchanges <-> r-i and r% «-> r%. 




3 Linear stability analysis 



By using the small-fc expansion method l|Rowlanddll969t llnfeld Il985t llnfeld and Rowlands! EfflOoh , we 
now investigate the linear stability of periodic waves to long-wavelength perturbations with wavevector 
fc(cos 9, sin 9 cos ip, sin 9 sin ip) where 9 is the angle between the direction of the wavevector and the x-axis, 
and ip is the azimuthal angle. We start from the ansatz 



u = uq(x) + e<f>(x)e ik ( x cos e+y sin 9 cos v+z sin 9 sin * 3 )~ ia ' < 



(3.1) 



where uq{x) is a periodic solution to (jl.ip . e«l, and the eigenfunction <&{x) must have the same period 
as Uq(x). Substituting (|3.1[) into (jl.ip and linearizing with respect to e gives 



da; 



L$ = iw$ - ifccos^Q* - 3ikcos9$ xx + k 2 (l + 2 cos 2 6>)$ x + ifc 3 cos6»$ 



(3.2) 



where 



da; 2 



— u + bu^ 2 — V, 



and $ and w are written as the expansions, 



$ = <J>o + k$i + 

uj = w\k + L^ik 2 



(3-3) 
(3.4) 



For the remainder of the calculation we follow a similar procedure to that first given in iParkesl (|l993l) . 
After substituting (|3 . 3[) and (|3.4p into (|3.2p and equating coefficients of k n we obtain the sequence of 
equations 

= i? M (a;) (3.5) 

in which the expressions for R nx are of the same form as in Part 1 after replacing jj by — iu)j. Since 
Luq x = 0, the solution to £$ n = R n + B n , where B n are integration constants obtained on integrating 
(US]), is 

= u 0x v n (3.6) 

where 

i r r x i 

(3.7) 



— |^4 n + f (Rn(x') + B n )uo x (x) dx' 



and A„ are additional constants. On integrating (|3.7p . secular (non-periodic) terms will occur in v n . To 
remove these we must insist that 

(v nx ) = (3.8) 

where 

(/}=Fpi ff(x)dx, (3.9) 



A is the period of uq, and Fp stands for Hadamard's finite part (jZemanian Il965h . Equation (13. 8|) provides 
a relation between A„ and B n which can later be used to help eliminate these constants. 
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To lowest order in k we have (L&o) x — 0. As a result of the translational invaria nce of up, this has a 
solutio n proportional to uq x . This result can be obtained more explicitly, as is done in lMunro and Parked 
(IHBSi), by using the consistency conditions (vq x ) = 0, ((L$i) x ) = 0, and (uq(L^i) x ) = to show that 
vox = 0. Without loss of generality we choose a unit constant of proportionality (which corresponds to 
setting vq = 1) and we are left with 

$ - u 0x . (3.10) 

Integrating the first-order equation gives 

L^i — iojiUQ — 2icos6 uq xx + Bi, (3-11) 

and after using (|3.7|) one obtains 

1 (a , iw i u o 



vix = —o- M + -^r 1 -icos6> Uo x +B lUo . (3.12) 



Then applying (|3.8p results in the relation 



2 

in which we have introduced the quantities 



MPq + - icos6» + Bi/3i = (3.13) 



'lie 



0* = \dr)- (3-14) 



After using (13. 6|) and (13. lip , the second-order equation may be written as 

(L<f> 2 ) x = iw 2 Wte + iuJiu 0x vi + u 0xx + ui cos9u - iB 1 cos6» - 2i cos 9 (u 0x Vi) xx . (3.15) 

To obtain u>\ it is not necessary to evaluate $2- Instead, we first apply the finite-part averaging operation 
(•) to (|3.15p . After using partial integration to show that (uo x vi) = — (uqvi x ), and by virtue of the 
periodicity of $2 (which implies that ((L<fr 2 ) x ) = 0), we get 

— koi (uqVi x ) + u>i cos 9 ai — i-Bi cos 9 = (3.16) 

where we have defined 

a s = (u s ) . (3.17) 

We then multiply (|3 . 1 5[) by uq and apply (•). The left-hand side can be shown to be zero by integrating 
by parts and then using the self-adjoint property of L and the fact that Luq x = 0. This leaves, after 
further manipulation via partial integration, 

(ugUix) — ( u ox) + w i cos9a 2 — \B\ cos#ai + icos6* (uq x Vi x ) = 0. (3.18) 

From (|3.12p we can obtain 

(uqVix) = MPi + 1UJ ]^ 3 - i cos 9 ai + Bi/3 2 , (3.19a) 

(ulv lx ) = + + (3.19b) 

(ul x v lx ) = A x + 12 -icos9(uo x ) + Biai, (3.19c) 
and after replacing u in (|2.2p by u and applying (•) we have 

/ 2 \ „ t r 86 1 

\ u 0x) = Ca l + Va ? ~ 15 a5 / 2 ~ 3 a3 - 

Substituting (|3. 19a.|) into (|3 . 16[) and (]3 . 18[) and then eliminating A\ and B\ from these two equations and 
(|3.13p leaves the following equation for u>\ 

a + a\U\ + a 2 uo\ + a 3 uu\ — 0, (3.20) 
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where 



a-o = (A) ( u lx) sin2 # + cos2 ^) cos ®i 
01 = sin 2 0, 



a 2 = ( /?i/?3 - - A &\ - 7/?o/?4 ) cos 



a 3 = i(/3 1 2 /34 + /3 2 3 + /3o/3 3 2 -/3o/32/34-2/3 1 /3 2 /3 3 ). 

Owing to the fact that uq x is zero at some points, the direct evaluation of the /3 S would require a 
finite-part calculation. This can be avoided by instead expressing these quantities in terms of the a s . To 
accomplish this, we require a number of recursion relations. The first of these is obtained by multiplying 
(12. ip by Uq/uq x , applying (•), and then simplifying the left-hand side using partial integration. This yields 

sa s -i = -p. + V(3 S+1 - y/3 s+3 /2 - \l3s+2- (3.21) 



Applying the same procedure to (|2.2|) and then replacing s by s — 1 gives 

15 A+3/2 ~ 3 



a s _! = Cf3 s + Vp s+1 - ^/3 s+3/2 - \Ps+2- (3.22) 



Eliminating /3 s _|_3/2 from the above two equations gives 

9C{3 S + 3V/3 S+1 + p s+2 = 3(5 - 4s)a s _i, (3.23) 
and putting the values s = 0, 1, 2 into (|3.23p generates the following three equations involving the required 

9C/3 + 3V/3i + p 2 = -15a_i, (3.24a) 
9Cf3 1 + 3Vf3 2 + = 3, (3.24b) 
9C/3 2 + 3Vf3 3 + (5 A = -ax. (3.24c) 

A further two equations for the f3 s are found by first eliminating /3 s + 2 from (|3.2ip and (|3.22[) to give 

30C/3, + 15VP S+1 - 4bf3 s+3/2 = 15(3 - 2s)a s _ x . (3.25) 

Putting s = 1 and s = 3/2 into this equation and eliminating f3 5 / 2 , and then using the resulting equation 
and (|3.25[) with s = to eliminate /3 3 / 2 gives 

156 fa + IbVCpx H — - = — Ca-i H — . (3.26) 

4 15 2 4 

Eliminating /3 7 / 2 and /?5/2 from the equations obtained from (|3.25ll with s = 1, 2, 5/2 yields 

UCA + uvca + - . H c _ 2tQs/2 _ 1^1. (3 . 27) 

4 15 2 4 

Using (|3.24a|) . (|3.26p and (|3.27p . all the j3 s for s = 0, . . . , 4 can then be expressed in terms of cti, a_i and 

«3/2- 

We now turn to the evaluation of a s . A recursion relation involving only a s can be obtained by 
multiplying (|2.ip by ug and (|2.2p by sug -1 , adding, and then averaging. If uq > or if s ^ 0, the average 
value of Uq will be finite and equal to a s and we may then write 











a- 


- |) « S+ 2 = 















2 8.s 

3 + 15 



ba s+3/2 . (3.28) 



However, if uo(x) is zero at some values of x, as is the case for the type I pc and II pc solutions, and s < 0, 
the average of ufj will no longer be finite, and in cases where the coefficient of an infinite integral is zero, 
(|3.28p has to be modified. Before continuing, it should be noted that, in contrast, (|3.2ip is always valid 
for s — since for this value of s the left-hand side originates from {uqxx/uq^ which is identically zero. 
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From (|3.28p it is evident that we will have to evaluate at least two of the a s directly. The simplest to 
find, due to the fact that the periodic wave solutions are of the form ug(x) — [r(x)] 2 , are oix/2 and a\. The 
evaluation of these integrals for type I, II pc , and I pc solutions is given in the Appendix. To determine the 
(3 S and (uq x ) we also require a_i, 0:3/2, «2, 0(5/2, an d ct3- Putting s = — 1 in (|3.28|) presents no problem 
as the only coefficient that is zero is multiplying a term originating from a finite integral. Thus we have 



1 / 4bai 



/2 



3C V 15 

To find a 3 / 2 we need to use s = —1/2. In this case (I3.28P must be re-written in the form 

3^«i/2 66ai 3(1 + 2s) 

0*3/2 = — + Jn }/2 Fp ^ u dx. (3.29) 

For type I solutions, the final term on the right-hand side is zero. For the type I po and II po solutions, 
the integral in (|3.29| is infinite and the finite part would have to be found numerically. In such cases 
a 3 / 2 needs to be obtained directly, as is done in the Appendix. The remaining a s can be found in a 
straightforward manner by putting s = 0, |, and 1 into (|3.28p which give, respectively, 

a 2 = C + 2V ai fl, 

ZCa 1/2 9Va 3/2 2lba 2 

0:5/2 = 



2 4 15 

9C* 12V 366a 5/2 
a3 = _ ai + _ a2 ___ 

The values of a s corresponding to finite integrals obtained using the procedure outlined were checked 
by numerical integration for specific values of the parameters. The numerical values of the remaining 
quantities, namely, a_i for the periodically extended solutions and the /3 S , for which the fin ite-part 



operation is not redun dant, were checked using a finite-part numerical integration technique (|Q'Keh 
1993t |Phibanchonll2006l ). 



4 Growth rate of instabilities 

Having obtained the three roots to the nonlinear dispersion relation (|3.20|) , we discard the real parts as 
they are of no importance in the context of stability. The solution is unstable if two of the roots are 
complex conjugates. If lo\ is one of these roots, the first-order growth rate of the instability is given by 

7 = 71 fc = I ImcJil k. 

When examining the dependence of the growth rate on the type of solution and the direction of the 
perturbation we find it convenient to introduce the parameter c, a rescaled version of C, defined by 

C (4.1) 



\Kr + )\ 



provided that V > — 66 2 /25. Then, if V > 0, the linear limit corresponds to c = — 1 and the soliton limit 
occurs at c = 0. The type I pc and II pc solutions have c > 0. As in Part 1, we only consider the stability 
of solutions for which b > as these are the more physically relevant. 

We feel that a plot of 71 against 8 shows the angular dependence of the growth rate more clearly than 
the more traditional approach of using a polar plot to depict the dependence of the real and imaginary 
parts of u> at all angles. The value of 71 for type I solution instabilities as a function of angle for a 
number of values of c and b are shown in Fig. [3] For the soliton limit (when c = 0) the growth rate is 
proportional to sin 9 which is in agreement with the results of Part 1. For the cnoidal wave solutions 
(when — 1 < c < 0), 71 is only non-zero above a critical angle, C rit, which increases with decreasing c. It 
is also evident that # max , the angle at which the maximum growth rate occurs, differs from 90° for cnoidal 
waves. The variation of both max and 9 CI1 t with c is shown in Fig. Sfa). The growth rate is largest for 
the soliton limit. From the plot of 7i m , the maximum value of 71 (the value when 8 = f? max ), in Fig.[5ja), 
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Figure 3: Plots of 71 against 9 for type I solutions with V = 4 and c of (top curves), —0.25, —0.5, 
-0.75, and -0.99 (lowest curves), (aj b = (b) b = 2 (c) b = 4 (d) b = 50. 




Figure 4: Plots of # max [solid lines], 9 CI a [dashed lines], and 9 C [dotted lines] against c: (a) V = 4 with 
b taking the values (top curves), 2, 4, and 50 (lowest curves); (b) V = — 4 with 6 taking the values 5 
(outermost curves), 10, and 20 (innermost curves). 




Figure 5: Plots of 7i m against c: (a) V = 4 with b taking the values (top curves), 2, 4, and 50 (lowest 
curves) - dots indicate value of 71 calculated for solitary wave solutions in Part 1; (b) V = — 4 with b 
taking the values 5 (top curves), 10, and 20 (lowest curves). 



it is apparent that there is a rapid variation in growth rate as c approaches zero. This is not unexpected 



8 



given that the waveform period increases rapidly and becomes infinite at the soliton limit, c = 0. Notice 
that the results found for the soliton limit are in agreement with the analytical results given in Part 1 of 
this study. 

In Part 1 it was found that 71 for solitary waves decreases with increasing b for a fixed value of 
r\. As is apparent from equation (2.2) of Part 1, for fixed rj, the amplitude decreases as b increases. 
However, if the amplitude is fixed (by using the appropriate value of rj in each case) then it is found 
that 71 increases with b. Cnoidal waves for different values of b but with the same amplitude and values 
of m will have different periods. It therefore seems inappropriate to compare the growth rates in such 
cases. Nevertheless, meaningful comparisons can be made on examining the entire growth rate curve as a 
function of 9. As can be seen from Figs. [3] and 03a), there is a more marked variation of the growth rate 
with 9 for angles above 9 max as b increases, and the plots in Fig. 0Ja) indicate that # ma x deviates from 
the perpendicular most of all when b is large. On the other hand, # cr it shows only a slight dependence 
on b. 

We now turn to the stability results for solutions in the form of functions extended by periodicity. 
When V > and c is increased above zero, one obtains type II pe solutions. It can be seen from Fig. [6](a)- 
(c) that the first-order growth rates of these solutions are higher than that for the soliton limit at some 
angles, but this range of angles decreases with increasing b. In addition to an increasing # cr i t with c, 
the first-order growth rate for exactly perpendicular perturbations vanishes for large enough b and c. 
Evidently the growth rate has a significantly greater angular dependence than for the Type I solutions. 




Figure 6: Plots of 71 against 9 for type II pe and I pe solutions. In (a)-(c), V = 4 and c takes the values 
(leftmost curves), 0.5, 1.0, 1.5 and 2.0 (rightmost curves) and (a) b = (b) b = 4 (c) b = 50. In (d) 
V = —4 and c takes the values 50 (top curves), 10, and 1 (lowest curves) with 6 = 5 [solid lines], b = 10 
[dashed lines], and b — 20 [dotted lines]. 

For the stability results we have examined so far, the first-order growth rate is non-zero for angles 
just below 90°. In the case of type I pe and II pe solutions when V < the results in Fig. EJd) indicate 
that there is a cut-off angle 9 C above which 71 vanishes. Hence such waves are, to first order, stable to 
perpendicular perturbations. In addition, the instability occurs over a relatively small range of angles, 
even for large values of c. For these types of solution, as is shown in Fig.[5fb), the growth rate increases 
monotonically with c, in contrast to the behaviour near the type I to type II pe transition. There is no 
spike in the growth rate at c = 1, the type I pc -II pe transition, since there is no sudden change of period 
around that point. 
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5 Conclusions 



This article has dealt with the small-fc stability with respect to transverse perturbations of cnoidal wave 
solutions to the SKdVZK equation which governs strongly magnetized plasma with slightly non-isothermal 
electrons. It was found that the growth rate of instabilities for ordinary cnoidal waves that travel faster 
than long- wavelength linear waves has a stronger angular variation as the distribution of electrons becomes 
increasingly flatter than the isothermal Maxwellian. We also examined a class of solutions that do not 
occur for equations without a square root term. These solutions, which are written as functions extended 
by periodicity have a minimum value of zero. This causes difficulties with some instances of recursion 
relations used in the determination of the growth rate, and results in it being necessary to evaluate an 
additional integral. This type of solution, for the case when the wave velocity is less than that of long- 
wavelength linear waves, has, to first order, a relatively narrow range of perturbation angles at which 
instability occurs and is stable to perpendicular perturbations. 

The half-integer nonlinear term in the SKdVZK equation was originally introduced by Schamel to 
model the effect of trapped particles in Bernstein-Greene-Kruskal (BGK) solutions of the Vlasov-Poisson 
equation. Hence Part 1 and this paper are to be viewed as a step towards the more formidable problem 
of studying the stability of the BGK modes where trapped particles play a significant role. 



Acknowledgements 

The authors wish to gratefully acknowledge support from the Thai Research Fund (PHD/47/2547). Two 
of the authors (S.P. and M.A.A.) also thank Warwick University for its hospitality during their visits. 



A The evaluation of 0:1/2, &i an d O3/2 

For cnoidal wa ye solutions, cti /?, Qa ■ and 0:3/3 are elliptic integrals and for their evaluation we therefore 
rely heavily on lBvrd and Friedman! (|1954l ) to which the result numbers in the following refer. 

Type I solutions 

The type I solution as given by (|2.5[) has a period of 2K(m)/r) where K(m) is the complete elliptic 
integral of the first kind. From result 340.01 we obtain 

, p\m) 

a 1/2 = r 1 + (n-r 1 ) ' p J A.l 
' K (m) 

in which II(n|m) is the complete elliptic integral of the third kind. 
Applying result 340.02 to yields 

_ (r-4-r-i) 2 f pE(m) + {p 2 + 2p(l + to) + 3m}U(-p\m) 
011 ~ 2(p+l)(l + m) ^ K{m) 9 ~ ™ 

|2|0 / n(-p|m) 

+ r 1 +2n(r 4 -n)-j^ r (A.2) 
where E(m) is the complete elliptic integral of the second kind. 

Type IIp C solutions 

After introducing 

A — B An- Br„ - /"* dX 

CTl 



A + B' An + Bru J (1 - acn(X\m)y 

the integrals ot n /2 for n = 1, 2, 3 may be written in the form 



12 x\ A- B J ^ \ p 



1 / Ar\ — Br u \ " v-^ f n 



" (a/a, - 1)%. 



10 



From results 341, 



h = ^(-g^lmj + ^tan-^sd^lm)] ) (A.8) 



where q = cr 2 /(l - cr 2 ), (j) = am(%|m), yu = y/m + q, sd(x|m) = sn(x|m)/ dn(x|m), 

- _ {2m - (2m - l)a 2 }h + a 2 E(4>\m) - {m + (1 - m)cr 2 }x + Ti 
2 ~ (l-<7 2 ){m+(l-m)<7 2 } ' 

- _ 3{2m - (2m - l)a 2 }I 2 + 2mF(4>\m) - {6m - (2m - l)cr 2 }/i + T 2 
3 ~ 2(1 -cr 2 ){fh+(l-fh)cT 2 } ' 

where F(^|m) is the elliptic integral of the first kind, and 

Y _ o- 3 sn(x|m) dn(x|m) 
™ (1 — (Tcn(x|m))" _1 ' 

Type Ipc solutions 

The a n / 2 for n = 1, 2, 3 may be written in the form 

n n / \ 

°V 2 = ( « ) ( r4 / ri _ x ) P/ p' 

where 

j = r ^ 

" 7 (l + psn 2 (X|m))™' 
From result 400.01, Ji = II(— p; (j)\m), where <fi = am(x|m). Results 336.01 and 336.02 yield 

_ gg(#n) - (m + p) x + {2p(l + ro) + 3m + p 2 }/i + Tfj 
2 ~ 2(1 + p)(m + p) 



(A.4) 



T _ mF{<j)\m) - 2{p(l + m) + 3m} Ji + 3{2p(l + m) + 3m + p 2 }I 2 + T 2 

4(l + p)(m + p) ' (A - 5) 

respectively, where 

= p 2 sn(x|m) cn(x|m) dn(x|m) 
(1 + psn 2 (x|m)) n 
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